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On the determination of the spin of the black hole in Cyg X-1 from 
X-ray reflection spectra 
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ABSTRACT 

The spin of Cygnus X-1 is measured by fitting reflection models to Suzaku data covering the 
energy band 0.9-400 keV. The inner radius of the accretion disc is found to lie within 2 grav- 
itational radii (rg — GM/c^) and a value for the dimensionless black hole spin is obtained of 
0.97^Q Q2^. This agrees with recent measurements using the continuum fitting method by Gou 

et al. and of the broad iron line by Duro et al. The disc inclination is measured at 23 -V^j'^ deg, 
which is consistent with the recent optical measurement of the binary system inclination by 
Orosz et al of 27 ± 0.8 deg. We pay special attention to the emissivity profile caused by irradi- 
ation of the inner disc by the hard power-law source. The X-ray observations and simulations 
show that the index q of that profile deviates from the commonly used, Newtonian, value of 3 
within 3rg, steepening considerably within 2rg, as expected in the strong gravity regime. 

Key words: accretion, accretion discs - black hole physics - line:profiles - X-rays:general 



1 INTRODUCTION 

Astrophysical black holes are characterised by just mass and spin. 
The measurement of spin requires observations which probe the 
immediate environment of the black hole event horizon. Two meth- 
ods which have emerged using cuiTent X-ray observations of lu- 
minous accreting black holes rely on fitting either the shape of 
the reflection component of the spectrum including a broad iron 
line, or the shape and the flux of the quasi-blackbody continuum. 
Both assume that the accretion disc extends in to the Innermost 
Stable Circular Orbit (ISCO), the radius of which is determined 
by the black hole spin. The Reflection Fitting method essentially 
measures the largest gravitational redshift of the disc, which comes 
from the ISCO, and yields that radius in units of gravitational radii 
{vg = GM/c^). The Continuum Fitting method measures the area 
of the inner disc and obtains the radius of the ISCO in km. Here we 



apply the reflection method to Suzaku X-ray data of the first stellar 
mass black hole, Cyg X- 1 . 

Early attempts to measure the spin parameter of Cyg X-1 us- 
ing the reflection method led to inconsistent results (Miller at al 
2005; Miller et al 2009). Adopting a Newtonian emissivity profile 
(surface flux on the disc varies with radius r as r^^ with ^ = 3) 
and using a special timing mode of XMM, a recent measurement 
by Duro et al. (2011) gives a dimensionless spin parameteiQ of 
a = 0.88^q '[[. Recent results using the Continuum Fitting method 
have been obtained by Gou et al (2011). They find a near ex- 
treme black hole with a > 0.95. The work relies on accurate op- 
tical determinations of the mass (14.8 ±0.1 M©) and inclination 
(27 ± 0.8 deg) obtained by Orosz et al (201 1). 

Here we study the reflection/iron line approach in detail, tak- 
ing account of the relativistic effects on the emissivity profile ex- 
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' The dimensionless spin parameter a = cj /GM^, where J and M are the 
angular momentum and mass of the black hole, respectively. 
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pected if the black hole spins rapidly. Such effects demand that 
the emissivity be steep if the disc extends to the innermost region 
within ~ 2rg (Wilkins & Fabian 2011). We find a = 0.97+jj [j^'^, 
which is a higher than, but marginally consistent with, that reported 
by Duro et al (201 1). Both reflection results are consistent with the 
Continuum Fitting result of Gou et al (201 1). 

We find the disc to be fairly highly ionized with an ionization 
parameter around 1800 erg cm s^', which means that the reflec- 
tion spectrum is dominated by a large edge due to ionized iron. 
This blurred edge carries information about spin in the spectral fit. 
Early work emphasised how this edge can be important in fitting 
the spectrum of Cyg X-1 (Ross, Fabian, Brandt 1996). 



2 SPECTRAL FITS 

The dataset we use is one of the 20 Suzaku observations of Cyg 
X-1 analyzed by Miller et al (2012; the observation of 2009 Apr 
8). The focus of that work is not on measuring the black hole spin, 
but rather studying the detailed disc-jet coupling in Cyg X-1. The 
spectra in that work have been fitted with a relativistically-blurred 
reflection model, but in order to make the work tractable, the fits 
kept the source inclination and iron abundance frozen at 27 deg 
and 1, respectively. The emissivity profile of the blurring function 
was also just a single power-law, fixed for the Newtonian value of 
3. We relax these assumptions here and apply them to an average 
dataset. We later tested the model on several other datasets (2009 
May 6, 25 and Dec 1) and obtained consistent results. 

The best fit to the 1-500 keV spectrum 
is shown in Fig. 1 and values of key parame- 
ters are given in Table 1. The model used, CON- 
stant *tbabs *mtable(windabs. fits) (gaussian-l-gaussian 
-i-diskpbb-i-powerlaw+kdblur2f*atable(extendx.mod)) 
HIGHECUT, is shown in Fig. 2. KDBLUR2F is a fast adaptation 
by JSS of the broken power-law emissivity convolver KDBLUR2. 
EXTENDX is a version of the self-consistent slab reflection model 
REFLIONX (Ross & Fabian 2005) extended to allow the photon 
index to be less than 1.4. (We do not use the REFHIDEN models of 
Ross & Fabian 2007, which incorporate a blackbody component 
because they are only available for Solar abundance iron and at 
the blackbody temperature found here for the low state of about 
0.16 keV the differences are minimal.) Absorption features due to 
the stellar wind of the companion (e.g., Hanke et al. 2009, Nowak 
et al., 2011) are modeled using WINDABS.FITS (Miller et al., 2012) 
and we refer to that paper for further details. The absoiption has 
little effect on our results. We note that the HIGHECUT model 
used to fit the data above 100 keV is phenomenological and is not 
strictly consistent with the 300 keV exponential cutoff assumed by 
the EXTENDX model, but has limited effect on the results presented 
here. The quality of the fit is very good, and the emissivity index 
q is found to be high > 6.8. The inclination of 23.7+j'^ deg is 
consistent with the optical measurement of Orosz et al (201 1). 

It is clear that the spectrum is reflection-dominated below 
10 keV (Fig. 2), although the whole source, from an energy point 
of view integrating over the whole spectrum, is not reflection dom- 
inated. This is a consequence of the hard iiTadiating continuum 
where much of the energy absorbed from Compton recoil due to 
photons at ~ 100 keV is emitted by the surface layers if the disc be- 
low 6 keV. The ratio of the 0.1-1000 keV flux from the reflection 
component is 0.70 times that from the power-law component, so 
the reflection fraction, calculated as the ratio of incident to emitted 
flux at the disc surface, !2 = 0.7. The blackbody disc component 
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Figure 1. Spectrum of Cyg X-1 in the Suzaku XIS0+XIS3 (black), PIN 
(red) and GSO (green) detectors, from Miller et al (2012). The ratio to the 
best fitting model (Table I) is shown below. The data have been rebinned 
for display purposes only. 
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Figure 2. Best-fitting spectral model plotted as EFe- The relativistically- 
blurred reflection component is shown as the red curve. The 100 keV kink 
introduced by HIGHECUT has no effect on our spin results. 



is 0.4 times the power-law and if reflection is responsible for some 
of that emission then ^ ~ 1. A narrow emission line is introduced 
at 6.4 keV to account for (weak) distant reflection (e.g. from the 
stellar companion) and a narrow absorption feature is included at 
1.85 keV to compensate for problems in the model response near 
the detector silicon edge. 

We show the model spectrum over the 2-65 keV band in the 
top panel of Fig. 3. The inner radius has been stepped from 1.3 to 
1.9, and then to 400rg and the innermost emissivity index q = 9. 
The effect of fixing ^ at 3 is shown in the lower panel: the model is 
now less sensitive to the innermost radii. 

We next fit the spectrum with a model in which the blurring is 
carried out on many small contiguous annuli, simultaneously. This 
is similar to the approach taken by Wilkins & Fabian (201 1) for the 
AGN IH 0707-495. Once again, the emissivity has a steep index at 
small radii before breaking to a flatter part, then dropping further 
close to the Newtonian expected value of 3. Simulations produced 
by ray tracing in the strong gravity regime revealed in Wilkins & 
Fabian (2012) show that this triple power-law shape is expected 
from an irradiated disc around a rapidly spinning black hole. The 
outer break suggests that the height of the irradiating source above 
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Table 1. Values of key model pai'ameters used in the spectral fitting. 



Component 


Parameter 


Value 


q fixed 


powerlaw 


Photon index, F 


1 ,7+0.014 
'■^'-0.011 


^■-'^-0.009 


kdblur 


Inclination, deg 

Index, qi 
Index, q2 


23 

1 60+°"' 
>6.8 

2.75±0.15 
4.0±1.1 


39.8l|-^ 
< 1.66 
3 
3 


reflionx 


Iron abundance / solar 
lonisation parameter, ^ 
erg cm s"' 


1.42l°;«4 
1765+?« 


0.89+0.12 
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Figure 3. (Top) Model reflection spectra, relativistically blurred for inner 
radii of 1.3, 1.5, 1.7, 1.9 and400 rg with an emissivity index (/of 9. (Bottom) 
Similar to above but with q = 3 and expanded scales. Note the reduced 
sensitivity to small values of rj^. 
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Figure 4. Emissivity profile obtained by fitting the blurred reflection model 
to small increments in radius. The profile has the characteristic steep - flat- 
moderate profile expected from irradiation by a source at a height ft = 5 — 
7rg above the black hole. The innermost annulus with a significant detection 
is at r ~ 1 .4r„ . 




Figure 5. Effect of varying the inner radius ri„ as a function of iron abun- 
dance Fe. 



the disc is around 5 — 7rg. The results show that emission from 
within 2rg is definitely required. 

The model which was used to fit the data from the XIS, 
PIN and GSO simultaneously uses KDBLUR2F, which is a broken 
power-law capturing the essence of Fig. 4. We investigated the de- 
pendence of the inner radius r[„ on various key parameters, the iron 
abundance and the inclination (Figs. 5~7). 

The data require an iron abundance which is slightly superso- 
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Figure 6. Effect of inclination on inner radius ri„ 
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lar, an inclination less than '--^ 30 deg and an intermediate ionization 
parameter around 1760. All results strongly point to an inner radius 
less than 2rg, even if q is fixed at 3 (Table 1, Fig. 7). 

As a test, we have fitted the PIN data alone, since this dataset 
is sensitive to the reflection hump. The iron abundance is fixed at 
unity (the depth of the edge in the XIS data indicates that the abun- 
dance cannot be either very low or very high). We also fix the in- 
clination at 27 deg, obtaining again < Ivg. This confirms that 
fitting the reflection hump alone can in principle measure spin. 

The reflection fits presented here all strongly point to the spin 
of the black hole being high. Using the broken power-law emis- 
sivity model model RELCONVF (a faster modification by JSS of 
RELCONV of Dauser 2010, 2012) yields a = 0.97+° °2'* at the 90 
per cent confidence level. The results for the other parameters are 
similar to those shown in Table 1 . 

We have also used the Novikov & Thome (1973; see also 
Page & Thome 1974) emissivity profile which is appropriate if 
the corona is immediately above the disc. The emission then varies 
with radius in the manner expected from a radiatively-efficient thin 
disc. The resulting fit has = 2422/2422, significantly worse than 
the earlier best fit. High spin {a > 0.98) is again required and the 
inclination is higher (/ 39 deg). 

2.1 The Newtonian power-law emissivity profile, = 3 

Duro et al (201 1) measure the broad iron line of Cygnus X-1 as seen 
by XMM-Newton in the EPIC-pn modified timing mode. They re- 
port that the black hole is spinning close to maximal with value 
a = 0.881q ][, on adopting the Newtonian, single power-law emis- 
sivity profile with index fixed at 17 = 3. A fit with q free gave 
q = lO^g and a = —0.1 ±0.4. They rightly dismiss the statistically 
more probable, high q - low spin, solution as making no physical 
sense, since high q is only applicable at high spin. 

Fitting our spectra with q fixed at 3 leads to a much worse 
(although acceptable) fit with A;|^^ — +194 (Table 1). The inner ra- 
dius is inferred to be at less than 1.66rg, the inclination is above 
35 deg (which is inconsistent with the optical result) and the ion- 
ization parameter is higher at nearly 3000. If the inclination is fixed 
at 25 deg, then < 1.9rg. The conclusion that the black hole is 
rapidly spinning remains robust. 

The spectral residuals obtained from using a single fixed value 
of q and a broken power-law with variable indices are shown in 
Fig. 7. The large difference in is mostly due to the red wing of 
the line feature at 2.5-3.5 keV and (to a much lesser extent) the 
edge at 6.5 — 7.5 keV. The residuals for fixed q are at about 2 per 
cent and increase to 3 per cent or more if is increased to 3rg and 
beyond. This demonstrates where our information is coming from 
and shows that from a statistical point of view we can probe the 
strong gravity regime. 

A single power-law emissivity has only limited validity. It is 
clear from the fitted emissivity profile in Fig. 4 that a slope of 3 
is a fair approximation from r = 2f-g outward. It underestimates 
the profile within r = Ivg and is therefore a poorer probe of the 
innermost region (see Fig. 3). If ^ = 3 is used, then it will likely 
yield an upper limit to the inner radius and thus a lower limit on the 
spin. Doubly broken power-law profiles are required when fitting 
for the inner radius and/or spin when the spin is high (rj^ < 2rg and 
a > 0.94). 

Cyg X-1 is a good example where a steep inner profile is 
needed (Fig. 8). In contrast, the source GX339-4, which has an in- 
ner radius of ~ 2rg and thus moderate spin, requires only a mildly 
steep index q of around 4 (Fig. 7). (We use the XMM data of the 
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Figure 7. Ratio of (heavily binned) spectral datapoints to the best fit- 
ting model with a single power-law emissivity profile (index q fixed at 3: 
black points) and to the broken power-law model with variable indices (red 
points). 
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Figure 8. Vaiiation of inner emissivity index qi on inner radius rjn for Cyg 
X-1 (top) and GX339-4 (bottom). A brolcen power-law emissivity profile is 
used with the outer index, q2, fixed at 3, for both plots. 
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Figure 9. Model spectrum with two power-law components. The index of 
the incident spectrum for the reflection component equals that of the steeper 
power-law. 

low-state observation of GX339-4 from Reis et al 2008 for this fig- 
ure.) 

The higher state of Cyg X-1 during the Duro et al (2011) ob- 
servations, in which the black body component contributes to about 
3 keV, may complicate variable q measurements. 



3 MODIFICATIONS TO ALLOW FOR THE TIME 
VARIABILITY OF CYG X-1 

In this Section we consider what modifications can be made to the 
spectral model developed here to allow for the complex time vari- 
ability observed in Cyg X-1. The aim is not to build a complete 
model which accounts for everything that has been observed but to 
indicate the directions in which modifications are needed. 

Study of the rapid variability of Cyg X-1 requires high count 
rates and thus instruments with large collecting areas, such as EX- 
OSAT (Belloni & Hasinger 1990), Ginga (Miyamoto et al 1988, 
1989; Gierlihski et al 1997) and RXTE (Nowak et al 1999; Revnivt- 
sev, Gilfanov & Churazov 1999; Gilfanov, Churazov & Revnivtsev 
2000; Pottschmidt et al 2003). The power spectrum is flat up to 
about 0.02 Hz, above which it drops with frequency / as up to 
a few Hz where it steepens to f^^. Time lags are seen over the fre- 
quency range 0.1-30 Hz in which soft bands lead hard bands, with 
the longest time lags being 0.05 s. As a light crossing time, this 
corresponds to a distance of ~ lOOOrg, which is much larger than 
the radius of any of the regions inferred by our spectral fitting. The 
lag timescale must therefore be some form of propagation time. In 
the low state, such as observed here, the iron "line" is found to be 
variable up to a frequency of a few Hz (Revnivtsev, Gilfanov & 
Churazov 1999). This is too low a frequency to be compatible with 
any light crossing time effects expected in our spectral model. 

Several models have been made for the spectral variability and 
time lags. Some early models are reviewed by Poutanen (2001). 
A model by Lyubarskii (1997) in which fluctuations are generated 
by and propagate inward through the disc has been extended by 
Arevalo & Uttley (2007) to explain many of the features observed. 
Most recently the hard state lags have been modelled as propagat- 
ing oscillations in a hot inner flow in the central region of an accre- 
tion disc truncated at about 20rg (Kawabata & Mineshige 2010), 
prompted by results of an earlier Suzaku observation of Cyg X-1 
reported by Makishima et al (2008). In that work the corona was 
modelled as having two optical depths, which allow for softer and 



Table 2. Values of model parameters obtained from the best fit with two 
power-law components of spectral indices Fs and Fi, . The incident spectrum 
for the reflection has index Fj. 



Component 


Parameter 


Value 


powerlaw 
powerlaw 


Photon index, F, 
Photon index. F^ 


1.45f 
1.25f 




Rintg 


1.63 


reflionx 


Iron abundance / solar 
lonisation parameter, ^ 
erg cm s^' 


1.28 
1775 




X-/doi 


2219/2419 



harder continua. Reflection was dealt with by a phenomenological 
model involving cold unblurred reflection and a broad gaussian iron 
line. 

The results obstained from the fits presented here reveal an 
inner radius of less than 2rg strongly irradiated by a corona situ- 
ated at ---^ 5 — 7 rg . This is incompatible with an extensive inner 
corona. The single power-law continuum used here does not allow 
for time lags unless it varies in slope on the required timescales. 
To allow for slope variations we have fitted the spectrum with two 
power-law continua, the steeper of which gives rise to the reflec- 
tion. The spectral fitting is degenerate provided that the sum of 
the two power-laws approximates the single power-law shown in 
Fig. 2. For illustration purposes we have selected a pair of spectral 
indices which give a good fit (Fig. 9; Table 2). Reflection continues 
to dominate the spectrum below 10 keV. 

A simple physical interpretation arises if the source of the 
steeper power-law is quasi-static and situated close to the centre of 
the disc, therefore producing the reflection. The harder component 
could originate from further out within the moving jet; both soft 
source photons from the disc and reflection produced by the con- 
tinuum are reduced by beaming (Beloborodov 1999). Time lags 
then occur if variations propagate up the jet from the inner, qua- 
sistatic and softer component to the outer, jetted harder component. 
They need not propagate at jet speed, but could be due to changes 
in the magnetic structure or collimation of of the jet. The power 
in these components and the jet is presumed to derive from the 
accretion disc and be transferred by magnetic fields coupled to a 
range of radii, so varying in a manner which is consistent with the 
Lyubarskii model. 

The lack of rapid reflection variations reported by Revnivtsev 
et al 1999) could in part be due to light bending effects implicit 
in the proximity of the inner, softer, power-law component to the 
black hole (Miniutti & Fabian 2004), if some of the rapid variability 
is due to (small) changes in height of that inner component. 

Further development of such models is beyond the scope of 
this paper, which is focussed on the low state X-ray spectrum of 
Cyg X-1. 



4 DISCUSSION 

We have used the X-ray Reflection Fitting Method to obtain a ro- 
bust measurement of the dimensionless spin parameter of the black 
hole in Cyg X-1 of a = 0-97^0 02^. Our result agrees with the con- 
tinuum fitting value of Gou et al (2011) and yields an even higher 
spin than the reflection measurement of Duro et al (201 1). Cygnus 
X-1 definitely appears to have high spin. 



© 0000 RAS, MNRAS 000, 000-000 



6 A. C. Fabian et al 



10 



9 



7 



6 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Spin 



Figure 10. Emissivity - spin contours for (left, black) a simple power-law 
profile and (right, red) a broken power-law profile with the outer index fixed 
at 3. The minimum x~ for the fit using the broken power-law is about 150 
below that for the fit using a single power-law emissivity profile. The results 
in Table 1 were obtained with the outer index free. 



We argue that non-Newtonian values must be used for the in- 
ner emissivity profile of the reflected emission of the innermost disc 
when the spin is high (rjn < Ivg and a > 0.94). 

The Reflection Fitting method rehes on detection of the in- 
nermost radius of the disc which is dense enough for reflection to 
occur. We assume that this is the ISCO, for within that radius the 
plunge orbits mean that the gas density drops rapidly so the reflec- 
tion signal vanishes. Although magnetic fields may have some ef- 
fect here, we consider it reasonable to assume that the dense parts 
of the disc that give the reflection spectrum are not seriously af- 
fected (see discussion and simulations in Reynolds & Fabian 2008 
and Shafee et al 2008). There is nevertheless a systematic uncer- 
tainty here which simulations need to resolve. 

The reflection fraction ^ ~ 0.7 is low compared with the sim- 
ple expectation of a source close to the black hole at /i ~ 5rg. Light 
bending should give ,^ ~ 2 (e.g. Fabian et al 20 11). The observed 
radio emission (Miller et al 2012) means that a jet is operating, 
so it is plausible that the emission model is more complex, as dis- 
cussed in Section 3. We suggest there that irradiation of the disc 
by the quasi-static base of the jet at a 5 — Ivg gravitational radii is 
responsible for the reflected emission, whereas the observed power- 
law continuum is dominated by higher, faster parts of the jet which 
mildly beam emission in our direction. 

Systematic uncertainties in the models and (relative) spectral 
calibration can also affect the details of our fits and thus results, par- 
ticularly of the innermost region. The Continuum Fitting method is 
also subject to systematic uncertainties in both models and spectral 
calibration (relative and absolute). 

Wide agreement on the high spin of Cygnus X-1 boosts con- 
fidence in spin measurements from the Reflection Fitting method 
of other objects as well as AGN, where it is the only method yet 
available. 

The spectral fits strongly indicate that the innermost radius of 
the disc lies within 2rg. Relativistic effects are strong there and a 
simple Newtonian assumption for the emissivity profile is inappro- 
priate. The effects are inescapable and involve both special relativ- 
ity, since the velocity of the reflecting matter on the disc is high, and 



general relativity, through gravitational redshift and light bending 
(Wilkins & Fabian 2010, 2012). 

Theoretical emissivity profiles due to point sources at differ- 
ent heights on the rotation axis above the disc plane are shown in 
Fig. 12 of Fabian et al (201 1). When the source height exceeds 3rg, 
the emissivity profile breaks at radius r ^^hio the Newtonian value 
of 3 outside and flatter within. The inner profile rapidly steepens to 
a high value within 2rg, if the disc extends in that far. 

We recommend that ^ = 3 be used initially. Should this lead 
to high spin being suspected then a power-law broken at least once, 
with a high iimer index (^i > 3), should be used. This is illustrated 
in Fig. 10 where we have fitted our Cyg X- 1 data by a single power- 
law emissivity of free index. It shows the best fit being of high index 
and relatively low spin {q ~ 10, a --^ 0.6), which is rejected by the 
physical sense argument of Duro et al (2011). If instead we use a 
broken power-law emissivity fixing the outer index q2 = 2>, then we 
find that the inner emissivity q\ remains high and the spin rises to be 
very close to maximal (a > 0.95). This now makes physical sense. 
Allowing the outer index to be free gives the results of Table 1. 
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